Hot-electron effect in spin dephasing in n-type GaAs quantum wells 
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We perform a study of the effect of the high in-plane electric field on the spin precession and 
spin dephasing due to the D'yakonov-Perel' mechanism in n-type GaAs (100) quantum wells by 
constructing and numerically solving the kinetic Bloch equations. We self-consistently include all 
of the scattering such as electron-phonon, electron-non-magnetic impurity as well as the electron- 
electron Coulomb scattering in our theory and systematically investigate how the spin precession 
and spin dephasing are affected by the high electric field under various conditions. The hot-electron 
distribution functions and the spin correlations are calculated rigorously in our theory, ft is found 
that the D'yakonov-Perel' term in the electric field provides a non- vanishing effective magnetic field 
that alters the spin precession period. Moreover, spin dephasing is markedly affected by the electric 
field. The important contribution of the electron-electron scattering to the spin dephasing is also 
discussed. 

PACS numbers: 72.25.Rb, 72.20.Ht, 71.10.-w, 67.57.Lm, 73.61.Ey 
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I. INTRODUCTION 

Spintronics is an active field which studies processes 
that manipulate the spin degree-of-freedom of electrons, 
with the goal of developing new electronic devices with 
improved performance and new functionality compared 
to the traditional ones which are based on precise control 
of the charge distribution of electrons. ^'^ Understanding 
spin dephasing is an important prerequisite for the real- 
ization of such devices. As most of the semiconductor 
electronic devices are very small and even a small ap- 
plied voltage gives a strong electric field, these devices 
usually work in the hot-electron condition.'^'"' Therefore 
understanding spin dephasing in the presence of a strong 
electric field is of particular importance to the spintronic 
application. 

Recent experiments have shown that the elec- 
tron spin lifetime is very long in rt-type Zinc-blende 
semiconductors. ^'^'^ In theory the spin dephasing in 
semiconductors without high electric field has been ex- 
tensively studied. Three spin dephasing mechanisms 
have been proposed:* the Ellit-Yaffet mechanism^do 
which is important in the narrow-bandgap and/or high 
impurity-doped semiconductors; the Bir-Aronov-Pikus 
mechanism^ ^ which is important in the pure or p- 
type semiconductors; and the D'yakonov-Perel' (DP) 
mechanism^^ which is the main spin-dephasing mecha- 
nism in n-type Zinc-blende semiconductors such as GaAs 
and InAs. The DP mechanism originates from the spin- 
orbit interaction in crystal without the inversion center 
and results in spin splitting of the conduction band at 

^ 0. This is equivalent to an effective magnetic field 
acting on the spin, with its magnitude and orientation 
depending on the electron wavevector. Moreover, an im- 
portant many-body spin dephasing mechanism due to 
combined effects of the inhomogeneous broadening of the 
spin precession and the spin conserving scattering by ir- 



reversibly disrupting the phases between spin dipoles has 
been proposed recently^^ and closely studied. ^"'■^^'^^'^'^ 

The study of the effect of electric fields on electron 
spins in semiconductors just begins. Experiments have 
shown that the spin polarization is not destroyed by 
the strong applied electric field in transport up to a 
few kV/cm.^*d9 jg revealed that under right config- 
urations the electric field can drive the electrons to ap- 
proach a larger spin injection length. ^■^'^'^^'^^•^'^'^^'^^'^^ In 
Ref. 27, the spin dephasing in quantum wires under high 
electric field is studied through Monte-Carlo simulation. 
The electric manipulation of the spin of two-dimensional 
(2D) electrons through the Rashba^* spin-orbital inter- 
actions using the in-plane AC electric field has also been 
proposed. Nonetheless how the hot electron effect af- 
fects the spin dephasing/transport is so far not fully in- 
vestigated. A complete understanding of the hot-electron 
effect on the spin dephasing in n-type GaAs quantum 
wells (QW's) can be obtained by solving the many-body 
kinetic Bloch equations^^^^'*'^^ which have been applied 
successfully to studying the spin dephasing^^ and spin 
transport recently. '^-'^ 

In this paper, we use the many-body kinetic equations 
to study the effect of the high electric field on the spin 
dephasing. The paper is organized as following: In Sec. 
II we present the model and construct the kinetic Bloch 
equations. Then we show the effect of the electric field 
on the spin dephasing problem by numerically solving 
the kinetic equations. In Sec. III(A) we first discuss how 
the electric field affects the spin precession. Then we de- 
vote ourselves to the understanding of the effect of high 
electric field on spin dephasing under various conditions, 
such as at different impurity densities; temperatures; and 
initial spin polarizations. We summarize the main results 
in Sec. IV. In Appendix A we present the effect of the 
electron-electron Coulomb scattering on the spin dephas- 
ing. 
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II. MODEL AND KINETIC EQUATIONS 

We start our investigation of an n-typc (100) GaAs 
QW of width a with its growth direction along the z- 
axis. An uniform electric field E and a moderate mag- 
netic field B are applied along the a;-axis (Vogit configu- 
ration) . Due to the confinement of the QW, the momen- 
tum along the z-axis of electrons is quantized. There- 
fore the electrons are characterized by a subband index n 
and a two-dimensional momentum k = [kx, ky), together 
with a spin index a{= ±1/2). For simplicity, we only 
consider the wells of a small width so that the separa- 
tion of the subband energy is large enough and therefore 
only the lowest subband is populated and the transition 
to the upper subbands is unimportant. It is noted that 
due to the so called "runaway" effect, the single subband 
model is valid only when the electric field is less than a 
few kV/cm.^^'^'* This is because when the electric field is 
above the threshold value, electrons gain energy from the 
field faster than they can dissipate it by emitting phonons 
and therefore the transition to upper subbands becomes 
significant. Consequently in the present paper we only 
study the case with the electric field up to 1 kV/cm which 
is sufficiently large to produce the hot-electron effect. 

For n-type samples, spin dephasing mainly comes from 
the DP mechanism. ^'^^ With the DP term included, the 
Hamiltonian of the electrons in the QW is given by: 

H = J2 {(£k-eE-R)5,,, 

+ [5MbB + h(k)] • ^}4^cka' + Hi. (1) 

Here = k^/2m* is the energy spectrum of the electron 

with momentum k and effective mass m* . cr are the Pauli 
matrices. R = {x, y) is the position. h(k) represents the 
DP term which serves as an effective magnetic field with 
its magnitude and direction depending on k. It is com- 
posed of the Dresselhaus term'^^ and the Rashba term.^* 
For GaAs QW, the leading term is the Dresselhaus term 
which can be written as: 

/ix(k) = ik,{kl - {kl)) ; 

hy{k) = ^kyiikl) - kl) ; 

/iz(k) = . (2) 

Here {kl) represents the average of the opera- 
tor — {-§^Y o'^®'^ electronic state of the low- 
est subband and is therefore {ir/a)^. 7 = 
i4/3){m* /mcv){l/ ^2m*^Eg){r]/ ^1 - 77/3) and rj = 
^/{Eg -h A), in which Eg denotes the band gap; A rep- 
resents the spin-orbit splitting of the valence band; and 
rricv is a constant close in magnitude to free electron mass 
rriQ.^^ The Rashba term is proportional to the total elec- 
tric field. For narrow band-gap semiconductors such as 
In As, the Rashba term is dominant. Whereas for wide 
band-gap semiconductors like GaAs, it is marginal in the 



regime of the applied electric field we study. The interac- 
tion Hamiltonian Hj is composed of the Coulomb inter- 
action Hf.e, the clectron-phonon scattering Hph^ as well 
as the electron-impurity scattering Hi. Their expressions 
can be found in textbooks. ^'"'^^ 

In order to study the hot-electron effect on spin de- 
phasing, we limit our system to a spacial homogeneous 
one in order to avoid the additional complicity such as 
charge/spin diffusion. The kinetic Bloch equations in 
such a system are constructed using the nonequilibrium 
Green function method with the gradient expansion^^ 
and can be written as: 

Pk.o-CT' - eE • VkPk.CTCT' = Pk.CTCT' Icoh + Pk.CTCT' Iscatt , (3) 

where pko-o-' represent the single particle density matrix 

elements. The diagonal terms describe the electron dis- 
tribution functions pk.o-cr = /ko-- The off-diagonal ele- 
ments i_ 1 = pI _ii = Pk stand for the inter-spin- 

'22 '22 

band polarizations (spin coherence).'^" The second terms 
in the kinetic equations describe the momentum and en- 
ergy input from the electric field E. p^^a' |coh on the right 
hand side of the equations describe the coherent spin pre- 
cession around the applied magnetic field B, the effective 
magnetic field h(k) from the DP term as well as the effec- 
tive magnetic field from the electron-electron interaction 
in the Hartree-Fock approximation. PkCTo-' Iscatt denote 
the electron-impurity, the electron-phonons, as well as 
the electron-electron scattering. The expressions of these 
terms are given in Appendix B. 

The initial conditions at t = are taken to be Pk(0) = 
and electron distribution functions are chosen to be 
those of the steady state under the electric field but with- 
out the magnetic field and the DP term. Specifically 
/k,CT(0) is the steady solution of the kinetic equations (3) 
with the spin coherence pk) the magnetic field and the 
DP term set to be zero. This initial distribution fimc- 
tions can be approached by assuming that at time —to 
there is no spin coherence Pk(— ^o) = and the electron 
distributions are just the Fermi distribution functions for 
each spin a at the background temperature T: 

hai-to) = {exp[(£k - Ma)/T]+l}"' . (4) 

and then self-consistently solving the kinetic equations 
(3) with the magnetic field and the DP term turned off 
(therefore no spin precession and pk = 0). By taking to 
to be large enough one may get the steady state solution 
before t = 0. In Appendix A we present a typical electron 
distribution function in the steady state under the elec- 
tric field. The imbalance of the the chemical potential 
Ml/2 7^ 1^-1/2 gives the initial spin polarization: 

^ iVe 1/2(0) -A^e -1/2(0) 

1/2(0) + A^e -1/2(0)' ^' 

where N,,„{t) = J^k fk,a{t) is the number of the electrons 
with spin a at time t. 
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III. NUMERICAL RESULTS 

As one notices, all the unknowns to be solved appear in 
the coherent and the scattering terms nonhnearly. There- 
fore the kinetic Bloch equations have to be solved self- 
consistently to obtain the electron distribution and the 
spin coherence. 

We numerically solve the kinetic Bloch equations in 
such a self-consistent fashion to obtain the temporal evo- 
lution of the electron distribution functions fuait) and 
the spin coherence pk(^)- Once these quantities are ob- 
tained, we are able to deduce all the quantities such as 
electron mobility n and hot-electron temperature Te for 
small spin polarization P as well as the spin dephasing 
rate for any spin polarization P. The mobility is given 
by A* = X)kcr/k<T(0) k/-^e) the electron temperature is 
obtained by fitting the Boltzmann tail of the electron 
distribution function; whereas the spin dephasing rate is 
determined by the slope of the envelope of the incoher- 
ently summed spin coherence p{t) = |Pk(OI-^^''^°'^^ It 
is noted that the spin dephasing time obtained in this 
way includes both the single-particle and the many-body 
spin dephasing contributions. 
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TABLE I: Material parameters used in the numerical calcu- 
lations 



We include the electron-electron, the electron-phonon 
and the electron-impurity scattering throughout our 
computation. For electron-phonon scattering, as we con- 
centrate on the relatively high temperature regime, only 
electron-longitudinal optical (LO) phonon scattering is 
important. The numerical scheme of the solution of the 
kinetic equations is laid out in detail in Appendix B. The 
total electron density Ne, the width of the QW a and 
the applied magnetic field are taken to be 4 x 10^^ cm~^, 
15 nm and 4 T respectively. The material parameters 
of GaAs are listed in Table I. The numerical results are 
presented in Figs. 1 to 5. 

A. Electric field dependence of the spin precession 
frequency 

In Fig. 1 we plot a typical temporal evolution of the 
electron densities of spin up and down for a GaAs QW 
with initial spin polarization P = 2.5 % under two elec- 
tric fields E = 0.5 kV/cm and E = -0.5 kV/cm at 
T = 120 K. S = 4 T for both cases. The corresponding 
incoherently summed spin coherence is also plotted in the 
figure. One can see from the figure that, the temporal 
evolutions of the electron densities and the spin coherence 
are similar to those in the absence of the applied electric 




FIG. 1: Electron densities of spin up and down and tlic inco- 
herently summed spin coherence p versus time t for a GaAs 
QW with initial spin polarization P = 2.5 % under different 
electric fields E = 0.5 kV/cm (solid) and E = -0.5 kV/cm 
(dashed). Top panel: S = T; Bottom panel: B = 4 T. 
Note the scale of the spin coherence is on the right side of the 
figure and the scale of the top panel is different from that of 
the bottom one. 



field. The electron densities and the spin coherence os- 
cillate as electron spins undergo the Larmor precession 
around the total (effective) magnetic field. Due to the 
spin dephasing, the amplitude of the oscillation decays 
exponentially. An interesting effect of the high in-plane 
electric field on the spin precession is that there is marked 
difference in the precession frequency under differcnit elec- 
tric fields (even the electric fields of the same magnitude 
but in the opposite directions). As shown in Fig. 1, al- 
though there is almost no difference in the corresponding 
spin dephasing rates, the periods of the oscillations are 
51.2 ps and 33.6 ps for applied electric field E = —0.5 
kV/cm and 0.5 kV/cm respectively. Both periods devi- 
ate from 40.6 ps, which is the electric-field-free one of the 
Larmor precession under the magnetic field B = 4 T. 

Moreover, it is expected that at very low temperature 
[i.e., a few Kelvin) where the momentum collision rate 
is small, the DP term can result in a rapidly damped 
oscillations in the spin signal when B = 0. At higher 
temperatures, due to the higher collision rates these oscil- 
lations disappear totally and the spin polarization decays 
exponentially with time^° and the oscillations can only 
be seen when there is an applied magnetic filed in Vogit 
configuration. Nevertheless, it is of particidar interest to 
note in the top panel of Fig. 1 that even at temperature 
as high as 120 K, the spin signal oscillates with period 
219.9 ps when there is no applied magnetic field but an 
applied electric field E = 0.5 kV/cm. 
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Both features above originate from the apphed high 
electric field E along the x-axis. With the applied electric 
field, the electrons get a net center-of-mass drift velocity 
Vd and the distribution function is no longer first-order- 
momentum free, i.e., X^k'^/ko- — m*Vd ^ 0. From Eq. 
(2) one finds that there is a net effective magnetic field 
B* along X-axis from the DP term which does not exist 
when E = ^. Prom the period of the spin oscillation in 
Fig. 1, one can deduce the effective magnetic field B* . 
When electric field E = ±0.5 kV/cm, the net effective 
magnetic field B* = =f0.74 T. 



# 

pq 




FIG. 2: The net effective magnetic field B* from the DP 
term versus the drift velocity Vd at T = 120 K with impurity 
density Ni = 0. The solid curve is the corresponding result 
from Eq. (7). 

The average of the total effective magnetic field the 
electrons experience at low spin polarization can be given 
approximately by: 



B -I- B* = B -I- 



1 /rfk(Ai/2-/k-i/2)h(k) 
gUB /rfk(/ki/2 - /k-1/2) 



(6) 

By taking the electron distribution function to be the 
drifted Fermi function /kcr = {exp[(k — TO*Vd)^/(2m*) — 
Mcr)/ (fcsTe)] + the effective magnetic field for small 
spin polarization can be roughly estimated as following: 

B* = 7m* ^V4Ef/[2{l - e-^f/"-^')] - Ec}/gf^B , (7) 

with Ef and Ec standing for the Fermi energy and con- 
finement energy of the QW respectively. In Fig. 2 the 
effective magnetic field B* deduced from the frequencies 
of our numerical result is plotted as a function of the 
drift velocity Vd for the impurity free sample. The re- 
sult predicated by Eq. (7) is also plotted in the same 
figure for comparison, with the hot-electron temperature 
Te obtained by fitting the Boltzmann tail of the calcu- 



lated electron distribution functions. It is seen from the 
figure that Eq. (7) gives a reasonable estimation of B* . 



B. Electric field dependence of tlie spin depliasing 
time of electrons with small spin polarization 
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FIG. 3: The SDT Tb versus the applied electric field E at (a) 
T = 120 K and (b) T = 200 K for initial spin polarization 
P = 2.5 % with different impurity densities: •, Ni = 0; 
Ni = 0.1 iVe; and A, Ni = 0.2 N^. 



In addition to affecting the spin precession frequency, 
the applied high electric field also changes the spin de- 
phasing time (SDT), i.e. the inverse of the spin dephas- 
ing rate, although the electric field does not couple to 
the electron spin directly. In Fig. 3(a) and (b), we plot 
the SDT of the electrons with initial spin polarization 
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P = 2.5 % as a function of the electric field E for dif- 
ferent impurity densities iV, at T = 120 K and T = 200 
K respectively. It is seen from Fig. 3(a) that for the 
impurity-free sample, the SDT first increases with the 
electric field from Ts = 175 ps at S = and then satu- 
rates to Ts = 300 ps, 70 % higher when E approaches to 
1 kV/cm. For samples with impurities, the SDT also in- 
creases with the electric field but with decreased increase 
rates for higher impurity densities. When the impurity 
density rises to 0.2 Ng, first increases for small elec- 
tric field and then decreases when the electric field is 
higher than 0.75 kV/cm. Moreover the change of the 
SDT with the electric field is much smaller than that 
of the impurity-free sample. The electric field depen- 
dence of the SDT varies as the background temperature 
changes. When the temperature is raised to a relatively 
high one, say T = 200 K in Fig. 3(b), the SDT increases 
slightly with the electric field and then decreases when 
E > 0.25 kV/cm even for the impurity-free sample. 

The electric field dependence of the SDT is understood 
due to the concurrent effects of the high electric field and 
the DP term. The most obvious effect of the electric 
field is that the electrons get a center-of-mass drift ve- 
locity and the center of the distribution functions drift 
away from k = 0. One consequence of the drift is that 
the DP term gives a net effective magnetic filed as dis- 
cussed above. This field is moderate and hence has little 
effect on the SDT. Another one is that because more 
electrons are distributed at large momentum region, the 
contribution from the DP term with large momentum is 
enhanced and the SDT can be reduced. Nevertheless, 
in addition to the drift, the high electric field also has 
another counter effect: As the high electric field gives 
the hot-electron effect with the electron temperature Tg 
higher than T, the scattering is strengthened. This can 
enhance the SDT.^'^^ In short, the drift of the center-of- 
mass in the momentum space trends to reduce the SDT 
while the hot-electron effect helps to enhance it in the 
regime of our study. With these two effects considered, 
the electric field dependence of the SDT can be under- 
stood. 

When the electric field is small, its effect on the DP 

term due to the drift is marginal. Therefore the SDT 
increases with the electric field due to the hot-electron 
effect when the temperature T is relatively low. As the 
electric field increases, the effect of the drift becomes im- 
portant and the SDT saturates consequently. It is noted 
that the hot-electron effect is more pronounced for the 
system with smaller impurity density under a given elec- 
tric field. As a result, the SDT increases slower with 
the electric field when the impurity density is higher. 
For high impurity-doped samples, the hot-electron effect 
is markedly smaller than that of the pure ones, there- 
fore the SDT only increases slightly in the low electric 
field region and then decreases as the effect of the drift 
dominates. Moreover, when the lattice temperature T 
increases, the hot-electron effect is also reduced. There- 
fore, in high temperature regime the drift effect becomes 



important even for low electric fields and it becomes pos- 
sible that the SDT may drop with the increase of the 
electric field even for the impurity-free QW's. Moreover, 
the change in Tg-E curve at high temperatures is smaller 
than that at low temperatures as shown in Fig. 3. 

It is noted that the electric-field dependence of Tg we 
obtain is different from that of quantum wires where the 
SDT decreases with the electric field. This difference 
may come from the different contributions of drift and 
hot-electron effect in quantum wells and quantum wires. 
In quantum wires as the electrons are much easier to 
be accelerated by the electric field to higher momentum 
states. Therefore the drift effect is more pronounced and 
it is possible that the SDT is reduced by the electric 
field. While the competing effect of the drift and the 
hot-electron in the QW results in a more complicated 
dependence of the SDT on the electric field. 
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FIG. 4: SDT Ta versus the electron temperature Te at lattice 
temperature T = 120 K with initial spin polarization P = 2.5 
% for Ni=0 (•) and iV, = 0.1 (♦). The curves are plotted 
for the aid of the eyes. 



In order to further elucidate the effect of the high elec- 
tric field to the SDT, we replot the the SDT as a function 
of the electron temperature Tg with T = 120 K for iVj = 
and Ni = 0.1 Ne in Fig. 4. It is seen that the SDT in- 
creases with the the electron temperature Te, similar to 
the electric-field-free case where the SDT increases with 
the temperature. ^^'^^ The figure also shows that the im- 
purities reduce the hot-electron effect and increase the 
SDT. These results indicate that for the electric fields 
we study, the electric field dependence of the SDT is af- 
fected mainly by the hot-electron effect in stead of the 
drift effect. 
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C. Electric field dependence of the spin dephasing 
time of electrons with high spin polarization 
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FIG. 5: The spin dephasing time Ts (•) vs. the initial spin po- 
larization for T = 120 K and A''^ = (a) and A^i = 0.1 A^e (b) 
under different electric field: E = 0; o, E = 0.25 Kv/cm; 
♦ , £ = 0.5 Kv/cm; A, E = 0.75 Kv/cm. 

We now turn to study the effect of the electric field 
on the spin dephasing with high initial spin polarization. 
Similar problem in the absence of electric fields has been 
studied in our previous work.^^ 

The numerical results are presented in Fig. 5 where 
the SDT is plotted as a function of the initial spin po- 
larization under different electric fields. It is seen from 
the figure that for all the electric fields we study, the 
SDT increases with the spin polarization as the case of 
E = 0.^'' However the speed drops with the increase of 



the electric field. It is interesting to see that the electric 
field dependence of the SDT is quite different for differ- 
ent initial spin polarizations. In low polarization regime, 
the SDT increases with the electric field while in high 
polarization one, it decreases with the electric field. For 
moderate spin polarized electrons, the SDT is insensitive 
to the electric field. 

The rise of the SDT with the initial spin polarization 
is understood due to the effective magnetic field from 
the HF term. As one component of this effective field is 
along the 2;-axis, it removes the "detuning" of the spin flip 
between the spin-up and -down bands and consequently 
suppresses the spin precession around the magnetic field 
and greatly reduces the spin dephasing. Therefore, all 
the factors, such as magnetic field, temperature, impu- 
rity, electron density, and applied electric field which can 
change the HF term, affect the spin dephasing in the high 
spin polarization case dramatically. These factors except 
the one of the electric field have been discussed in detail 
in our previous work.'^^ As for the factor of the applied 
high electric field, both the drift and the hot-electron ef- 
fects affect the HF term. The drift of the center of the 
mass in momentum space provides two competing effects 
on the HF term: One is to enhance the HF term through 
the net effective magnetic field B* discussed above. The 
other is to destroy the HF term by increasing the DP 
effect. Meanwhile the hot-electron effect tends to soften 
the HF term through the increase of the electron temper- 
ature and the scattering rate. Our results indicate that 
the electric field tends to reduce the effective magnetic 
field from the HF term in high spin polarization regime 
and consequently reduce the SDT. 



IV. CONCLUSION 

In conclusion, we have performed a systematic investi- 
gation of the spin dephasing due to the DP mechanism 
in the present of high electric fields by constructing a 
set of kinetic Bloch equations for n-type semiconduc- 
tor QW's based on the non-equilibrium Green function 
method with gradient expansion. In our theory, we in- 
clude the in-plane electric field, the magnetic field in the 
Vogit configuration, the DP spin-orbital coupling and all 
the spin conserving scattering such as electron-phonon, 
electron-non-magnetic impurity as well as the electron- 
electron scattering. By numerically solving the kinetic 
equations, we study the evolution of electron distribu- 
tion functions and the spin coherence of spin polarized 
electrons. The SDT is calculated from the slope of the 
incoherently summed spin coherence. In this way, we are 
able to study in detail how the spin precession and the 
spin dephasing are affected by the electric field in various 
conditions, such as the impurity, temperature, and spin 
polarization. 

The in-plane electric field has two competing effects on 

electron spins. The most obvious one is that the electrons 
get a center-of-mass drift velocity and the center of the 
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distribution functions drifts away from k = 0. One con- 
sequence of the drift is that the DP term contributes a 
non-vanishing net effective magnetic field which changes 
the period of the spin precession. The larger the electric 
filed is, the larger the drift velocity and consequently the 
net effective magnetic field is. For the electric fields we 
study, the net effective magnetic field is up to the order 
of 1 T. This moderate magnetic field has marginal ef- 
fect on the SDT although it results in a distinct change 
in the spin precession period. Another consequence of 
the drift is that because more electrons are distributed 
at large momentum regime, the contribution of the DP 
term with large momentum is enhanced. Therefore the 
drift can reduce the SDT. In additional to the drift, the 
high electric field also introduces another counter effect 
on the spin dephasing: The scattering, which tends to 
drive the electrons to the steady state, is enhanced as the 
hot-electron effect brought by the high electric field with 
the electron temperature T^. higher than the background 
one T. That is, the high electric field can also affect the 
spin dephasing through hot-electron effect. With these 
two effects of the electric filed on spin dephasing, the elec- 
tric field dependence of the spin dephasing is very rich in 
detail. 

In small spin polarized regime, the hot-electron effect 
tends to enhance the SDT as the increase of the scatter- 
ing rate reduces the inhomogeneous broadening. There- 
fore in the small electric field regime where the effect of 
the drift is marginal, the SDT increases with the electric 
field due to the hot-electron effect. For larger electric 
field, the effect of the drift become stronger. Therefore 
the SDT saturates under the joint actions of the drift and 
the hot-electron effect. When the impurity density or the 
background temperature T increases, the hot-electron ef- 
fect reduces and the effect of the drift becomes relatively 
important. As a result, the increase of the SDT with 
the electric field is reduced. For some regimes, the SDT 
decreases with the increase of the electric field when the 
drift effect dominates. 

In the high spin polarized regime where the HF term 
plays an important role in the spin dephasing, the hot- 
electron effect tends to reduce the SDT as both the in- 
crease of the electron temperature Tg and the increase of 
scattering reduce the HF term. Therefore, in the high 
spin polarization regime, the SDT decreases with the in- 
crease of the electric field. 
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FIG. 6: A typical distribution function of spin-up electrons 
in the steady state with electric field E = —0.75 kV/cm and 
P = 2.5 %. 



APPENDIX A: EFFECT OF COULOMB 
SCATTERING ON SPIN DEPHASING 

It is stressed that the electron-electron Coulomb scat- 
tering is of particular significance in this investigation. 
It is not only because the Coulomb scattering is crucial 

in the build-up of the hot-clcctron temperature and the 
hot-electron distribution functions, but also because it 
has strong contribution to the spin dephasing with or 
without the electric field. With the Coulomb scattering, 
the electron distribution functions become smoother in 
the momentum space and electrons distribute more uni- 
form around the drift center as shown in Fig. 6. 

It is noted that the treatment of the electron-electron 
scattering in the present paper takes account of the full 
effect of the Coulomb scattering which is different from 
our previous work^'' where the Coulomb scattering is 
evaluated by replacing the distribution functions and the 
spin coherence in the scattering with the corresponding 
isotropic averages along the angle. In this way we are able 
to compute the electron distribution function more accu- 
rate under high electric field and to have the hot-electron 
effect included in our calculation. It is noted that in the 
absence of the applied electric field, the approximation 
we used before greatly reduces the CPU time and gives 
good qualitative results. However, this approximation is 
not good if one tries to get the results quantitatively. 

We further show the importance of the Coulomb scat- 
tering in the spin dephasing by plotting the SDT as a 
function of the applied electric field with and without 
the Coulomb scattering in Fig. 7. As shown in the figure, 
for electrons with small spin polarization, the Coulomb 
scattering tends to drive the electrons to the equilibrium 
state when E = Q or the steady state when applied with 
a finite electric field, and hence greatly reduces the in- 
homogeneous broadening originated from the DP term. 
As a result the Coulomb scattering increases the SDT 
with/without applied electric fields. For the high spin po- 
larized system, as discussed before, the effective magnetic 
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field along the z-axis from the HF term plays an crucial 
role in the spin dephasing and the Coulomb scattering 
tends to reduce this effective magnetic field. Therefore 
the SDT becomes smaller when the Coulomb scattering 
is included. 
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APPENDIX B: NUMERICAL SCHEME OF THE 
KINETIC BLOCH EQUATIONS 



In this appendix we describe the scheme of the numer- 
ical solution of the Bloch equations (3). We first rewrite 
the Bloch equations as following: 



/k,cr = eE • Vk/k,<T + /k.crlcoh + /k,cr|scatt , (Bl) 

Pk = eE • VkPk + Pklcoii + Pk| scatt • 

(B2) 



FIG. 7: Spin dephasing time in GaAs QW with (filled cir- 
cle) and without (open circle) the Coulomb scattering in- 
cluded under different spin polarization: (a) P=2.5 % and 
(b) P=97.5 %. 



The coherent terms are 



5/i 



k,CT 



dt 



coh 



-2a{[giJ,BB + hx{k)]lmpi^ + hy{k)Rep]^} + Aalin ^ VqP^+qPk, 



(B3) 



dt 



+^Z]'^q[(/k+qi - /k+q-i)Pk - Pk+q(/ki " /k-i)] • 



(B4) 



In these equations Vq = J2q^ Ko f q^+g^+K^ ] [-^(^g^)!^ '^i*^ '*o standing for the static dielectric constant and = 
6TTNee^/{aEf) denoting the screening constant. The form factor \I{iqg)\'^ = ir"^ siv? y / [y"^ [y"^ — tt^)^] with y = qza/2. 
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The scattering terms are 

-2TvNi ^ C/q(5(£k - £k-q) [fkai'i- " fk-qa) " ^{PkPk-q)] " 27r ^ Vq^(5(£k-q - £k + £k' - £k'-q) 
q qkV 



dt 



\ -27r glxS{Su - ek-q - ^qq,\) [Nqg^xifka " /k-q<T) + /k^(l " /k-qa) " Re(pkPk-q) 



(1 - /k-q<7)/ka(l - /kv)/k'-qa' + 2^kPk_q(/kV - fk'-qa') + -^Pk' Pu' -q{h-cia - ha) 

— {k k — q, k' <-> k' — q}, 



(B5) 



5pk 
5t 



scatt 



{'^ 49.A^(ek - £k-q - f^qq.A) [pk-q(/kl + /k-l) + (/k-ql + /k-q-i " 2)pk - 2N^q^x{pk - Pk-q)] 

qg^A 



+7riV, ^ C/25(£k - £k-q) [(/kl + /k-i)Pk-q " (2 - /k_ql - /, 



k-q 



~X]'^^q'^('^k-q-ek+£k' -ek'-q)( (/k-qiPk + Pk-q/k-i)(/k'i "/k'-qi +/k'-i " /k'- 
qk' ^ 

+Pk[(l - /k'i)/k'-q| + (1 - /k'-|)/k'-q-| - 2Re(p;;,pk'-q)] - Pk-q[/k'|(l - /k'_q|) 

+(1 - /k'-i)/k'-q-i - 2Re(pj;,pk'-q)]') |-{k ^ k - q,k' ^ k' - q} , 



(B6) 




FIG. 8: The discrete momentum space. 



in which {k k — q, k' k' — q} stands for the same 
terms as in the previous {} but with the interchanges k ^ 
k — q and k' ^ k' — q. In these equations gq.q^.x are the 
matrix elements of electron-phonon couphng for mode 
A. For LO phonons, g"^ j^q 
9f)]}|/(%)pwitha 



{47rar!j^/^/[V27I(g2 



eV/^/(2^L0)('^~ -'^0^)- i^oc is 
the optical dielectric constant and ^^lq LO-phonon 
frequency. Nqq^\ = l/[exp{Qqq^x/kBT) — 1] is the Bose 



distribution function of phonon with mode A at temper- 
ature T. = E,. {47r^ieV[«o((/' + g^)]} V(ig,)|2 is 
the electron-impurity interaction matrix element with Zi 
standing for the charge number of the impurity. Zi is 
assumed to be 1 throughout our calculation. 

For numerical calculation, one first turns the Bloch 
equations into discrete ones. To facilitate the evaluation 
of the energy conservation, i.e. the (5- function in the 
scattering terms, we divide the truncated 2D momentum 
space into N x M control regions, each with equal energy 
and angle intervals as shown in Fig. 8. The k-grid points 
arc chosen to be the center of each control region and 
therefore are written as: 



k„,TO = V2m*^(cos^m,sin^„ 



(B7) 



with En ^ {n + 1/2) AE, Om = mAO. Here n = 
0, 1, • • • , iV- 1 and m = 0, 1, • • • , M- 1 with En-i = -Ecut , 
the truncation energy, and 6m-i = {M — 1)2-k/M. 

In order to carry out the integration of the ^-function 
in the scattering term, AE has to be chosen to sat- 
isfy VLq = nLoAE [or Oq = ("lo + \)^E]. Under this 
scheme, the coherent terms and the scattering terms of 
the electron-impurity and the electron-phonon scattering 
can be divided into discrete ones directly. Nevertheless 
/kff and pk in the Coulomb scattering terms are not all 
on the grid points we choose. We approximate them to 
be the interpolation of the nearest grid points with the 
same energy. 

The driving terms should be treated with caution as 

the equations arc stable only for some finite differential 
schemes, such as forward differencing and central differ- 
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encing schemes. In this study, we use the forward dif- 
ferencing scheme. However, the usual expression of this 
scheme is based on the Taylor series expansion and is dif- 
ficult to apply to the polar coordinate system which we 
use in this work. This difficulty can be circumvented by 
the so called discrete conservation principle:*^ 



eE- Vk/k,<Tlk= 



^ (fk eE • Vk/k,a 



m*AEAe 



— — / ds eE • n fk <7 

m*AEA9 7aa„^ -""'^ 

> ^ / eE • n fk cr 

m*AEAe Z^Jf, , , 

\ \ ' Tji ~ ri'm' n'm' r 



m*AEA0 



(B8) 



Here f2„,m and dVtn^m are the control region which con- 
tains the grid point \<in,m and the corresponding bound- 



ary. In the last step of the above equation, the integration 
of the boundary is replaced by the summation over the 
first order quadrature on the four (or three if the control 
region is the neighbor of k = 0) sides of the boundary 



drin^m with 



and nJJ ™ standing for the length and 



the outward normal to the boundary Qnm H Qn'm'- In 
order to satisfy the request of the numerical stability, 
^n,m is chosen to be kn,m if -eE • nJJJJ^ > and lin',m' 
otherwise. 



It is noted that this choice of k makes our approach 

identify to the forward differencing scheme. The time 
evolution is computed by the fourth-order Rungc-Kutta 
method.'*'* The computation is carried out in a parallel 
manner in the "Beowulf" cluster. For a typical calcula- 
tion, it takes about 7.5 hours to get one SDT with 16- 
node AMD Athlon XP2800-I- CPU's when both N and 
M are chosen to be 32. 
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